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' Abstract 

o ■ 

' A scalar model of the glueball is offered. The model is based on the nonperturbative calculation 

, of 2 and 4-points Green's functions. Approximately they can be expressed via a scalar field. On the 

O ' basis of the SU(3) Yang-Mills Lagrangian an effective Lagrangian for the scalar field is derived. The 

' corresponding field equations are solved for the spherically symmetric case. The obtained solution is 

. interpreted as a bubble of the SU(3) quantized gauge field. 

(N : 

^ ! 1 Introduction 

>\ . . 

^\ , The nonlinearity of quantum chromodynamics leads probably to such objects as a hypothesized flux tube 

OO ■ filled with a longitudinal color electric field and stretched between quark and antiquark and a glueball 

(N ■ which is a blob of gluonic fields. There is a difference between flux tube and glueball: the first object is 

^ ' created by the quark and antiquark (the sources of the color field) the second one has not any sources - it 

\ is a blob of selfinteracting fields. The flux tube has the directed electro-color field but for the glueball it 

' is not clear: whether is it something like monopole or no. These objects are impossible in linear theories, 
for instance in Maxwell electrodynamics, as its existence is connected with a nonlinear fields interaction 

^ • among themselves. 

I \ At the moment there are different glueball models such as flux tube [Jj , bags [2] , constituent gluons P] , 
QCD Hamiltonian in Coulomb gauge 0], non-Abelian Born-Infeld theory or the conjectured duality 

r~| ! between supergravity and large-N gauge theories 

' Here we present the glueball model where the gluon field is completely quantum one and it is described 

I by the nonperturbative manner. In this model the gauge field is organized by such a way that 

\ (^^) = t)ut {^(A^)^^ 7^ 0. In this case the color field in glueball performs nonlinear oscillations and 

' the nonlinearity of the Yang-Mills equations do not allow us to present these oscillations as quanta, i.e. 

the glueball is not a cloud of quanta. The basis for the presented glueball model is: (a) the initial SU(3) 
Lagrangian for the quantized field A^ is averaged over some quantum state \Q)\ (b) the 2 and 4-points 



Green's functions l^Af^A'^) and {A^A^A^A^^j arising in this case are approximated with help of 

some multiplet of scalar fields; (c) varying with respect to these scalar fields give rise to equations which 
describe these Green's functions. 



2 (J) approximation 

In any quantum field theory the Green's functions give us the full information about quantized fields. 
In this section we would like to present equations which will describe 2 and 4-points Green's functions 
by some approximate manner in the QCD. For this we will average the SU(3) Lagrangian where we use 
some approximate expressions for 2 and 4-points Green's functions. The SU(3) Lagrangian is 

£sui3) = \f^.F^'"' (1) 

where = - d^Aj^ + gf^^^A^A^ is the field strength operator; B,C,D = 1, . . . , 8 are the 

SU(3) color indices; g is the coupling constant; f^^^ arc the structure constants for the SU(3) gauge 
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group; is the gauge potential operator. In order to derive equations describing the quantized field we 
average the Lagrangian over a quantum state | Q) 

ig/^^^ - 9.1^ (x)) A^^{x)A^- [x)) + (2) 

\_^g2^BC.B,^BC.B. (A%^[^x)A^^^{x)A^^nx)A^-^{x)) . 

One can see that schematically we have the following 2, 3 and 4-points Green's functions: l(dAf' 



{(dA) A^) and {{Af)- At first we suppose that the odd Green's functions can be written as the following 
product 

(If (a;)A^(y)l^(^)) « (if (x)I^(y)) + (other permutations) = (3) 

as we have supposed that (Af (x)) — 0. Later we suppose that 2-point Green's function can be presented 
in so called one-function approximation pj as 

(If (x)A^(y)) = Ql'^ix^y) « (4) 

where <^^(x) is the scalar field which describes the 2-point Green's function. The 4-point Green's func- 
tion can be written in one-function approximation as the product of corresponding two 2-point Green's 
functions 



Af « (Af (A^(z)Af (u)) + 

Al{x)A^^{z)\ lA%{y)^{u)) + (if (^)) [A%{y)A^^{z)\ 



Taking into account these expression for the 2,3 and 4-points Green's functions we can derive an effective 
Lagrangian -Ce// — \ Cj for the scalar field 4>^ which describes 2 and 4-points Green's functions (for 
details, see Appendix El 

>Ce// = -\ {d,<i^^f + Y r r - rMf - Y m'^.f + 
-f [^^^"^ - cci' - f (CO' + irr) {r<t'n 

where the indices a — 1,2,3 are SU{2) indices, m — 4,5,6,7,8 are the coset SU{3)/SU{2) indices, (fi^ 
are some constants. The field equations are 

d^d^^r = ~r[2'^"''^'^ + ^i{rr-roro)], (7) 

a^a^^" = -0™ [2<^'^</.° + A2 (</.'"(/."- CO] ■ (8) 

One can note that in Ref. j7] was investigated a system of coupled scalar fields and was shown that such 
system may have soliton solutions by some specific choice of the potential term. 

3 Glueball as a bubble of 0^ field 

Now we would like to consider the spherically symmetric solution with the following ansatz for the scalar 
field 

r{r) = a = 1,2, 3, (9) 

(^'"(r) - m = 4,5,6,7,8. (10) 

V 10 



Let us note that this ansatz means that the components 0" have another behaviour then the components 
(j)™- . One can say that such situation is close to a colored fiux tube solution filled with the longitudinal 



electric field. After substitution ^ (|10|l into equation (jT)) (O we have 



+ = cl^[f + X,{ct>'~m')], (11) 

/" + -/ = f[<P' + X2{f-^^')] (12) 

where 20Q(/)g = and 2(f>^(j>^ = fp; m, fi are some constants which will be calculated by solving 
equations pil) and p2|) and we redefine Ai_2/2 Ai_2- Evidently these equations can not be calculated 
analytically. The preliminary numerical investigations show that this equations set do not have regular 
solutions by arbitrary choice of m, ^ parameters. We will solve equations Q) ||SJ) as a nonlinear eigenvalue 
problem for eigenstates 4>{x), f{x) and eigenvalues m, /i, i.e. we will calculate m, parameters such that 
the regular functions 4'{r) and /(r) do exist. 

At first we note that the solution depends on the following parameters: 0(0), /(O) and Ai^2- We can 
decrease the number of these parameters dividing equations H12|) on (/)'^(0). After this we introduce 
the dimensionless radius x = r(t){0) and redefine 4>{x) / (j){0) — > 0(a;), /(x)/0(O) f{x) and ra/(f){Q) — > m, 
n/4>(0) fi. Thus we have the following equations set 

+ V - <f>[f + X^{^'~m^)], (13) 

r + ^f - f['P' + X2{f-n')]. (14) 

We will search the regular solution with the following boundary conditions 

0(0) = 1, 0(oo) = m, (15) 
/(O) = /o, /((X3)=0. (16) 

One can say that (j>{x) is like to kink and f{x) to soliton. Let us rewrite equation 114|l in the following 
form 

-(/" + V)+/^e//=(A2M')/ (17) 

where we have introduced an effective potential 

K// = (0' + X2f) . (18) 

Immediately we see that with the boundary conditions (|16|l equation (|17|l is the Schrodinger equation 
and it may have a regular solution only if Ve/ / has a hole {(jp ^Zl^ const, -^ji^ q-j ^^^^ ^^^^ energy level 
A2/x^ have to be quantized. 



4 Numerical solution 

We choose the following numerical method for solving equations H14|) : we take a null approximation 
for the function f{x) (which is foix)) and solve equation (|13|l in the foUowng form 

0^' + ^^'o = 00 [fS + Ai - ml)] (19) 

where mo is the null approxiamtion for the parameter m, the boundary conditions are (|15|l and the 
function (j)o{x) is zero approximation for the function 4'{x). Having the regular solution 0o(a;) we can 
substitute it into equation H14|) and solve the equation 

/r + ^/( = /i[0^ + A2(/NA^?)] (20) 

with the boundary conditions (|16|l (for the numerical calculations presented here we take /q = -\/0.6). 
Thus we have the first approximation /i (x) which we substitute into equation p3|l 

0'/ + ^01 = 01 [/f + Ai (0? - m?)] . (21) 

This equation gives us the first approximation for the function 0i (x) and so on. On the i*'' step we will 
have 

0r + -0', = 0, [/^i + Ai (0? - m?)] (22) 



and 



(23) 



For every step we have the values mf and fif as an approximation for the true eigenvalues values m*^ 
and n*^. 

4.1 The more detailed description of the numerical calculations 

At first we will describe the numerical solution of equation (|19|l . For this we choose the null approximation 
for f{x) as 



fo{x) 



/U6 



cosh I 



(24) 



The typical solutions for the arbitrary value of mo are presented on Fig^ We see that by toq < mg (ttiq 
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Figure 1: The typical singular solutions for 
equation I^Hl- '^i =0.1. 



Figure 2: The typical singular solutions for 
equation ((201 • ^2 = 1-0. 



is an unknown parameter which gives us the regular solution) the solution (fioix) is singular and near to 
the singularity the equation has the form 

00 ~ (25) 

consequently the solution is 

M^) « \f^^— (26) 

where xq is some constant depending on mo. On the other hand by mo > ttiq the solution is presented 
on Fig^and the corresponding asymptotical equation is 



which has the following solution 



(f>'^{x) + -0;, w - (Aim^) 00 

sin (x^/Xim^ + a) 
(j)o[x) « 



(27) 



(28) 



where (f>oo and a are some constants. All of that allows us to assert that there is a value ttiq for which 
does exist an exceptional solution which with some accuracy is presented on Fig^] For this value m^ the 
equation H19|l has the following asymptotical behaviour 



cl)'^{x) + -<l>'o^2\, (mS)'(0o-mS) 



(29) 



and the corresponding asymptotical solution is 



(f>o{x) w mo + Po- 



(30) 



where 4>oo and (3 are some constants. 



The next step is finding the first approximation for the /i (x) function. The equation is 



(31) 



From the previous calculations one can assume that there is the exceptional regular solution 4'o{x) with 
the asymptotical behaviour (|30|l . Then the numerical investigation shows that for the arbitrary fi there 
are two different singular solutions which are presented on Fig|21 Analogously to equation (|25|l the 





Figure 4: The iterative functions </'i,2,3,4(a;). 



Figure 3: The iterative functions /i, 2, 3. 4(2^) 
singular behaviour of the function /i (x) is 

V M X - xo 

Evidently that we can suppose that there is a regular exceptional solution f*{x) by /ii 
following asymptotical behaviour 

g-x^ (m*)^-A2(Ml)^ 



A2 X - 

2" 



Xo 

1 



by fii < ^il, 
by fii > fil- 



(32) 
(33) 

/X* with the 
(34) 



where foe is some parameter. The next step is substituting the first approximation fi{x) into equation 
H13I) for finding the regular exceptional solution 4'l{x) by then (l)l{x) will be substituted into equation 
H14I) for finding the regular exceptional solution /|(a;) hy fi — 11 2 and so on. 

The result of these calculations is presented on Fig'sOl Eland Tabled We see that there is the con- 
vergence 0*(x) 4>*{x), f*{x) f*{x), m* m* and fi* fi* where f*{x), (ffix) are the eigenstates 
and m* , fi* are eigenvalues of nonlinear eigenvalue problem (|12|l . 

The verification of the presented numerical method was done for the soliton solution, for details see 
Appendix ^ 
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1.8374351... 


1.594328... 


1.6186108... 


1.61823766.. . 




1.492105312.. . 


1.4938287... 


1.4921473... 


1.4921473... 



Table 1: The iterative parameters m* and /i*. 



5 The properties of solution 

In this section we would like to describe the properties of the derived solution. It is easy to see that the 
asymptotical behaviour of the regular solution is 

-£c^2Ai(m»)^ 

cj)* (x) « m* + /? (35) 

X 

/* (x) « foe (36) 



where m* and ^* are the parameters derived in the coarse of iterative sohition of equations H13() H14() . 
The energy density of the presented solution is 



e{r) = — 
.9 



(37) 



here Ai_2,"T'* and ^* are redefined according the remark after eq. H12I) and we add the constant term 
— A2^^/4 for the finiteness of the full energy. Thus the glueball energy is 



oo 



47r 



00-^1 (Al,2,TO*,^* 



(38) 

here we have redefined r, / and </) according to remark before eq. (|13(l . The quantity (/)g ^ defines the 
radius of flux tube since the dimcnsionless variables x for flux tube is a; = p<j)^. The profile of the energy 
density is presented on Fig. The numerical calculations for the dimcnsionless integral Ii gives 



oo 



2 *2 

I — m 



2 A 



da; « 6.28. 



(39) 



Figure 5: The profile of the energy density. 



The next what we can do is the calculation the ratio of the energy glueball and string tension. On 
the basis of similar ideas presented here a flux tube solution was found in Ref. |S|. The basic idea 
is that the quantized SU(3) gauge potential A^,B = 1,2,..., 8 can be splitted on two pieces: (a) 
the potential components A^,a = 1,2,3 belong to the SU(2) subgroup, (b) the potential components 
A™, m = 4, 5, 6, 7, 8 are coset componets. For the flux tube case it is supposed that are almost classical 
degrees of freedom but componets (analogously to glueball case) are completely quantum degrees 
of fredom. 2 and 4-points Green's functions for the A™ components are expressed via a scalar field 
(what is similar to presented here glueball case) . The assumptions which are similar to glueball case lead 
to a numerical solution describing a flux tube filled with the longitudinal color electric field. The linear 
energy density (or string tension) is 



oc 



r 







dp - 



(40) 



— 0o/2 (A, 0oo,"li 2) 

5 



where 



A\{p)^lf; Alip)^^M, 03(^).M. 



(41) 



The numerical calculations give I2 ~ 0.63. 



The flux of the electric field is 



oo oo 

$ = y" Ef^ds = 27tJ pli^MPldp = y y xfix)v{x)dx ^ (A,0oo,mt_2) 



(42) 



where ~ fv/g is the longitudinal color electric field. Eq. I|42(l shows that (like to Coulomb law) the 
flux of the electric field is proportional to a color charge defined as q — 1/g. But of course there is a 
dimensionless correction I3 coming from the nonlinearity of the theory. The numerical calculations give 
I3 w 0.79. 

Let us consider the ratio 

K= (43) 

here we take into account the value of dimensionless constant /iir ^ 10. It is necessary to note that for 
the ratio (|43|l we consider the case when the scalar field cj) has the same magnitudes at the center of fiux 
tube and glueball. It can be compared with the lattice calculations where this quantity is presented 
as W/y/^ « 3.64 for 0++ glueball. 

It is interesting to consider the dimensionless ratio 



w 4/1 



0.04 (44) 



which tell us that the ratio H43|) is proportional to the flux of electric fleld. 

Now we would like to calculate the angular momentum of the glueball in the offered model. The 
angular momentum operator is 

M = I \fx Ie"^ X iJ^j 1 dV (45) 



where = Fg^ = doAf - d^A^ + y/^^^^i^if is the the operator of the color electric field; = 
erjkF^^^ is the operator of the color magnetic field, and F^f. = djA^ ~ dkAf + gf^^^AfA^; i,j,k = 
1,2,3. Let us consider 



E^ X 



One can show that the expectation value 



(46) 



{m,) = (47) 

as in the considered case either dQ{- ■ ■) = Q oi I^A^ Af • • • ^ = in the consequence of the presence the 

factor rj^i, in the assumptions Q and It means that the spin of the presented glueball model is zero. 
It is interesting to note that there is an opinion .L4, that pure glueball can only be spin 0. 

From this consideration immediately we see that in this approach the glueball with nonzero spin 
probably can be derived using ansatz similar but with nonzero correlation between A^ and Ap. 

Another interesting possibility for the future investigations is the derivation of a mass spectrum. For 
this we see two ways: the first one is the search of excited states on the basis of ansatz Q for the 2-point 
Green's function; the second one is the search an another ansatz for the 2-point Green's function which 
gives the glueball with another mass. The preliminary investigations show that probably excited states 
can not be derived using the presented iterative method. Thus the derivation of glueball mass spectrum 
in the nonperturbative approach is the complicated problem and it is the goal of the future investigations. 

Finally, we would like to touch upon the connection between our solution and Derrick's Theorem ^l5] . 
This theorem tells us that in 3 spatial dimensions the scalar field theory with nonnegative potential do 
not have absolute stable solutions with finite energy which means that at the infinity the solution must 
tend to a global minimum where V{global minimum) = 0. But the potential for the interacting scalar 
fields of © has global and local minima. Our solution is in one of the local minima. The Derrick's 
Theorem tells us that if we add two constant terms ^('^o'^o)^ ^^"-^ to the potential in eq. 

10) we will have the following potential 



V (global minimum) = 0, 

V (local minimum) > 0. 



(48) 



In this case the presented solution (in the full agreement with the Derrick's Theorem) will have the 
infinite energy. If we add only one constant term ^(0o"0o")^ we will have the following potential 

V (global minimum) < 0, ^^^-^ 

V (local minimmn) = 0. 

In this case the energy of the glueball solution is finite but the stability of the solution have to be 
investigated in the future works. 

After the discussion of the properties of derived solution one can see that the presented model is the 
model of glueball with zero spin where some combination of fields 0° and (jfm push out each other. 



6 Physical discussion and conclusions 



In this letter we discuss the glueball solution presented on Fig's|3l El This scalar model of glueball is 
derived with the assumptions that: (a) 2 and 4-points Green's functions of the SU(3) gauge potential 
can be approximately expressed via a scalar field; (b) the scalar fields components with a small subgroup 
indices belonging to SU{2) £ SU{3) can have the different qualitative behaviour in some physical situa- 
tions (in flux tube and glueball) in comparison with the scalar components which indices belong to the 
coset SU{3)/SU{2). As the consequence we see that in this case exists a blob of the quantized SU(3) 
gauge field and the coset components push out the SU(2) components of the scalar field that is like to 
the Meissner effect in superconductivity. Such solution can be interpreted as the glueball in a medium. 
It follows from the fact that 2 and 4-points Green's functions of 0™ are nonzero at the infinity. 

Remarkably that similar situation exists in a flux tube solution obtained in Ref . |H] . There is only one 
essential difference between flux tube and glueball solutions: in the flux tube solution the gauge potential 
components belonging to the small subgroup SU(2) in the first approximation can be considered as 
classical degrees of freedom that allows to exist a longitudinal color electric field directed from quark to 
antiquark. 

Let us to underscore that in this interpretation the derived bubble of the quantized SU(2) components 
live in the sea of the quantized coset components. But we can present an another interpretation of this 
solution in which both components exist in vacuum. 

One can unite two assumptions about 2-point Green's function and the term breaking the gauge 



mvariance 



'A^ix)Af,{y)) « -rj^p [/^^'^/^^^0''(a:)^^(2/) + /S^™/^^" (</>™(x)0"(y) - C^o)] 
Let us remind that the indices d,e— 1,2,3 and m, n = 4, 5, 6, 7, 8. In this case 



A^ix)A^{x) 



E 



[r^f {4>''{x)f + V (/'""")2 (0"(x) - 0^)' 



(50) 



(51) 



{r-'f{^\x)-^l) 



(52) 



which describe the variance of nonlinear oscillations of the gauge field and they are nonzero inside of 
the bubble only. It means that the quantized field is concentrated in this region. By such a manner the 
correlation between A^ (x) and A^ (x) , B ^ C components are nonzero in the same region. Consequently 
one can say that in this approach the quantized field SU(3) fields is concentrated in the bubble and can 
be interpreted as glueball in the vacuum. 

The presented here approach to the QCD is similar to a field correlator method with one difference: 
in our approach there is dynamical equations for the Green's functions which are derived from the SU(3) 
Lagrangian. 

This approach to the Green's functions which can be approximately considered as scalar fields (or a 
condensate) may have interesting applications for gravity where scalar fields have various applications: 
infiation, boson stars, non-Abelian black holes and so on. Our approach allows us to speculate that the 
nonperturbative quantum effects can be very important in some gravitational phenomenon. 
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A The effective Lagrangian 

In order to derive equations describing the quantized field we average the Lagrangian over a quantum 
state I Q) 



Q 



c 



I?/""" ((a,^? - A'^'A"') + js'/""'"' f"^'"' (a^-a^-a'^'-a''', ' 



Schematically we have the following 2, 3 and 4-points Green's functions: (^{dA) ({dA) A^) and (^{A) 
At first we introduce the 2-point Green's function 

{A^ix)A''piy))^g^^ix,y) (54) 

The first term on the rhs of equation (|53|l is 

(5^1? (x)) [d^A^'^ix]) = d^^d^ (If (a;)) [A^'^iy)) | = V^^d.^d^G^^ {x, y)\ (55) 



For the simplicity we consider the case with a;o — yo- For this Green's function we use so called one- 
function approximation |5] 

g^,^{x,y) « -,M/^^^/^^^0^(x)0^(y) (56) 

where (p^ {x) is the scalar field which describes the 2-point Green's function. Physically this approximation 
means that quantum properties of the field A^ can be approximately described by a scalar field (f)^ {x), i.e. 
in this approximation the Lorentz index n is not very important. Taking into account this approximation 
we have 

{[d,A^) (aM^^)) = -v:/''^^!''^'' = -12 (a^</.^) {d^^^) (57) 

and 

{{d^A^) [d^A^f^)) - -3 (9^0^) (9^0^) , (58) 
Later we suppose that the odd Green's functions can be expressed as the sum of the following products 
(^A^{x)A'^{y)A^{z)) « (^A^ix)A'j^{y)) {a^{z)) + (other permutations) = (59) 

as {A^{x)) = 0. It gives us 

{{dj^ix)) A'^{x)A^{x)) = d,, (If (x)l^(y)l^(z))|^ = (60) 

and consequently in our approximation 

( (d,A^ - a.lf ) l^'^I^"^) = 0. (61) 

For the last quartic term on the rhs of equation H53() we assume the following approximation 

[A^ix)Af,iy)A^iz)Af iu)) « (if (a:)l^(j/)) (l^(z)lf (^.)) + 

[A^ix)A^{z)) {A'^{y)Afiu)) + (A^ix)A^{u)) (l^(y)l^(z)) . 

In fact it is the assumption that 4-point Greens function is the product of two 2-points Green's function. 
In this approximation the Ihs of H62|l is 



A^{x)A^ix)A"'^ix)A^-iuy^ 



where Ai.2;(Pi ^ Qi 2) some parameter depending on the values of the indices -Pi, 2, Qi,2 

{Ai, if aU indices Pi, 2, Qi,2 = 1, 2, 3, 
A2, if aU indices Pi, 2, Qi, 2 = 4, 5,6,7, 8, (64) 
1 , otherwise 

where Ai,2 are some parameters. Introducing this index Ai,2;(Pi 2,Qi 2) would like to say that the 
presented approximate quantization procedure is a little different for the scalar field components belonging 
to the small subgroup SU{2) e SU{?,) and the coset SU {?,) / SU {2) . In this case 

Al (/^^^V'^-^'V"/) {/"^"'f^^'W'P'^) + A2 [fBEl^nj-CE,n^,n^n^ f^j D E2P j RE2q ^ ^q^ (65) 

(other terms) 

The calculations show that 

jABC j-ADR l^j^BEiPi jyCEiQi^Pi^Qi^^ f^jyDE2P2 J-RE2Q2^P2^Q2^j ^ Q, (66) 

^ABC ^ADR^jBE,ajCE,b^a^b^jDE2CjRE2d^c^d^ ^ (gy) 

jABC jADR ^jBEim jCEm^m^n^ ^jDE2P jRE2q^p^q^ _ q ^gg-j 

Consequently 

Ai,2;(P,2,Q.2) (/^^^^V^^^^^'/'^^(:^)0^^(:^)) (/^^^^V«^^^^'/'^n^)0^M^)) V^.V'" = 0. (69) 
The similar calculations show that 

jABCjADR ^jBE.P, fDE^Q^^P^^^^^Q,^^^~^ ^jCE2P2 J«£2Q2 (^)^Q2 (3,)) = 

1\ (70) 



27 



jABCjADR^jBE,ajDE^b^a^b-^^jCE2CjRE2d^c^d-^ ^ (71) 
jABC jADR (^jBEim jDEm^m^n-^ l^jCE2P jRE2q^p^q'^ _ _ ^ |-y2) 



Consequently 



Al,2;(P,2,Q,2) = 

07 07 07 (T3) 

-Al {rrf + yA2 (0"v")' + - irr) {r^ri ■ 

Analogously 

Al,2;(P,2,Q..2) (r^^^V^^^'^^'^^Hx)'/'^^^)) = 

27 9 27 9 27 (74) 

-Al {rct^i^ + yA2 (^'"^")' + - {rc^i irr) • 

Finally for x ^ y = z = u the quartic term is 

fARBfACD (^A^A^A^^'A°v'^ = yAi (^0")^ + ^As (0">'")^ + 81 (0"^) ('/>"0") • (75) 
Therefore we have the following effective Lagrangian describing 2 and 4-points Green's functions 



yAi irrf + yA2 (0"v™)' + 81 irr) ('^"V") 



(76) 

If we redefine 0° — > 20°/ (3g) and A1.2 — > Ai,2/2 we will have the ordinary Lagrangian for the scalar field 



Y^e// = {d^cf,^) + ^ (0^'^)' + ^ (0'"<^™)' + (<^>'^) (0'"0™) . (77) 

Now it is necessary to do an essential remark. The SU(3) Lagrangian for the gauge group is very 
nonlinear: it has terms. It is well known |lUj that in Xcf)'^ theory the similar nonlinearity give rise 



to an additional term to potential term. One can suppose that the similar situation takes place in this 
situation, too. Here we suppose that the nonlinear terms like leads to the appearance of some term 
in the initial Lagrangian. For the simplicity we assume that the mass term will appear. Thus the final 
form of the effective Lagrangian is 




where (t>Q are some constants. In this situation the field equations for the approximate scalar description 
of the QCD are 

S^a'^c^" = [2</.-</.- + Ai , (79) 

^^a^^'" = 0™ [20>° + A2 (<^'>" ~ CO] ■ (80) 

In conclusion we have to note that this procedure for the approximate calculations of 2 and 4-points 
Green's functions should be some approximation for an exact procedure which obtains all Green's func- 
tions bu a nonperturbative manner. At first such procedure was offered by Heisenberg for the quantization 
of a nonlinear spinor field [T^j and later was applied for the QCD [H]. 



B The numerical calculations of the soliton 

For the validation of the presented method of solving the nonlinear equations H12|) we choose a soliton 
solution. The corresponding equation is 

^ = y" = yi^-y')- (^D 

The solution is 

We rewrite the (|81|l equation in the form of the Schrodinger equation 

-y" + yV,ff^-Xy (83) 

where Ve/f = —y"^ and A = 1. It shows us that the regular solution exists only for a discrete spectrum of 
"energy level" A. We will solve this equation by an iterative procedure. At first we have the equation 

-y'i + yi i-yl) = -Xm (84) 

for the first approximation yi{x) and where Ai is the first approximation for the A. For the numerical 
solution we choose the null approximation as 

The typical solution for the arbitrary values of the parameter Ai is presented on FigEI This picture 
shows us that there is a value A| for which the solution is exceptional one. One can find this exceptional 
solution choosing the appropriate value of the "energy level" A*. After which an exceptional solution 
y* (x) is substituted into equation for the second approximation 1/2 (x) 

-2/2 - 2/2 {ylf - -A22/2 (86) 

and so on. The result is presented on Table[21and Fig[7| One can see that A* 1 and y*{x) is convergent 
to y*{x). 
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